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The dynamics of the 2D Potts ferromagnet when quenched below the transition temperature is 
investigated in the case of discontinuous phase transition, which is interesting for understanding 
the non equilibrium dynamics of systems with many competing equivalent low temperature phases, 
that appears to be not much explored. After briefly reviewing some recent findings, we focus on 
the numerical study of quenches just below the transition temperature on square lattices. We show 
that, up to a certain time, metastable states can be observed for which energy stays constant above 
the equilibrium energy and the self-correlation function displays a fast decay. 



I. INTRODUCTION 



The way in which thermodynamical systems attain equilibrium when driven through a phase transition has been 
subject of much work [3,[2, 3, 4]. In fact, while equilibrium requires generally an infinite time to be reached, different 
dynamical regimes can be observed which present often peculiar properties (different kinds of dynamic scaling, critical 
slowing-down, spatial self-affinity, etc.) depending on the thermal history and the system. Dynamics of different phase 
transitions have been widely investigated, e.g from a disordered to an ordered phase, or between two ordered phases. 
Among the model systems, the Ising model displays both instances: in zero magnetic field it undergoes a continuous 
phase transition from the high temperature paramagnetic phase to the low temperature ferromagnetic phase; the 
transition is discontinuous when at constant low temperature an external magnetic field is reversed. Both problems 
have been widely investigated. In the first case, one of the main results is contained in the Allen-Cahn law for the 
ordered phase domain growth at low temperature, which predicts a power law decay of the energy ~ t~ 1 ' k , with 
k = 2,3 [H, 0. For what concerns the second case, the metastability associated with the field-driven first order 
phase transition has been studied with nucleation theory and its generalizations 0, E[ > which describe effectively the 
nucleation time and other features of the metastable states [3, IE Is • 

Not many results seem instead available for the ordering dynamics of systems possessing many equivalent ordered 
phases [H, HI]. On one hand, the coarsening at low temperature presents peculiar properties not well described by the 
Allen-Cahn theory. This is predicted on general bases by the Lifshitz's criterium Q, according to which a system 
possessing many equivalent ordered phases may equilibrate towards non-homogeneous phases when quenched below its 
critical point if their number is larger than the real space dimensionality. Thus, equilibrium phases at low temperature 
can be not uniquely defined and systems may display complex behaviours. On the other hand, and differently from 
what happens for the Ising case, there exists little theoretical study of the metastability near the temperature-driven 
transition of the 2D q- Potts model (which is first order for q > 4). A recent one [l(| predicts in particular that, for 
this model, metastability is a finite size effect. The study is based on a droplet expansion and, to our knowledge, has 
not yet been related with the underlying dynamics [llT ] (metastability in the Potts model has been investigated with 
different aims in [ijj], and for the field driven transition in [Hj]). 

The aim of the present paper is to resume some recent findings on the dynamics of the d = 2 Potts model [14| and to 
present some new results concerning metastability in the undercooled states just below the transition. Investigation of 
temperature driven metastability in this model seems particularly important for several reasons. The model presents 
interest in itself, being an effective model for systems with many possible broken symmetries, and mimics the dynamics 
of a large variety of real systems, such as gas adsorption froth [l5[, grain growth [16| . biological cell [13, HH. On 
the other hand, as mentioned before, the existence, or not, of metastable states associated to the temperature-driven 
first order transition of the Potts model is an old problem, firstly proposed by Binder [ill ], and which remains still 
open. 

The next section quickly recall the Potts model and presents a brief review of some recent finding on its dynamics 
on the square lattice in different temperature regions below the transition temperature. After introducing some 
elementary notions about metastable states in thermodynamical systems and a brief discussion on some recent findings 
concerning the metastable states of the Potts model, Sec. 3 points out the presence of metastable states for the model 
on square lattice simulations in terms of energy and self-correlation function. Section 4 presents a short summary 
and perspectives. 
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II. POTTS MODEL DYNAMICS BELOW THE TRANSITION TEMPERATURE 

The Potts model is described by a Hamiltonian of the kind [HI : 

«{»7} = £(i-w. a) 

(v) 

where r)i is the state of site i that can assume one out of q different values, usually identified by an integer: rji £ [1, q]. 
When q — 2 Eq. fl} reduces to the Ising Hamiltonian (within multiplicative and additive constants). Here the case of 
nearest neighbour interaction on square lattice will be considered. It presents a discontinuous equilibrium temperature 
driven phase transition at T q for q > 4 fl4j |. 

The non equilibrium dynamics of the Potts model presents several peculiar features. For example the determination 
of the domain growth exponents in the Allen-Cahn law at low temperature remained problematic for a long time 
[l5l whereas it was just an artifact of the residual energy [2^, Hl| (see below). From a general point of view, a 
recent survey on the different dynamical regimes displayed by the Potts model has been done by Ferrero and Cannas 
[22| . These regimes are observed in finite-size realizations of the system when they are cooled below the transition 
temperature T q (via Monte Carlo dynamics) and depend on the final temperature T that, according to the dynamical 
features, can be grouped in roughly three different ranges [22], delimited by the temperatures < T g < T* < T c < T q . 

• For T* < T < T c relaxation is dominated by simple coarsening, which follows the quench according the Allen- 
Cahn law [l| for the growth of the domain size, withe energy per site decaying as : e(t) ~ £ -1 / 2 . 

• For T g < T < T* the simple coarsening is interrupted at long time scales and the system is trapped with finite 
probability in some configurations of high symmetry with characterisitic lengths that grow with the system size 
[23 ], These states can be identified with those predicted by Lifshitz, from which the system escapes by activated 
dynamics. 

• For T < T g the system gets stuck in some kind of frozen or glassy states with well defined lifetime, which 
diverges with T — > as e c / T [12] • An interesting aspect [U, [23|, [24[ is that during relaxation the system obeys 
a generalized Allen-Cahn law wihich includes an additive constant eo: 

e(t) ~ e + a ■ r 1/2 . 

The value of the constant corresponds to the average energy of the glassy states at zero temperature [24| ■ The 
existence of freezing at zero temperature and the role of activated process at low temperature was first pointed 
out by Vinals and Grant. Domains are pinned by some local finite energy barrier [25l . |26( which cannot be 
overcome if T ~ 0. For low finite temperature these states present a lifetime independent of the system size 
[22| , which is determined by the predominance of the activation time with respect to an equilibration time [25| . 

In addition to the above temperature ranges it is seen that, in the range just below the transition T q , some metastable 
states with energy and structure similar to the disordered phase can be observed [22l [2^ ] . Investigation of the nature 
of these states is the subject of the next section. 

III. METASTABLE STATES 
A. Theoretical premises 

Existence of metastable states just below a discontinuous transition is predicted by phcnomcnological theories like 
Van der Waals equation for the liquid-gas transitions [3] and are a general feature of theories with mean-field or 
long range interaction [J] . As in the Van der Waals equation, metastable branches appear as a continuation of the 
stable branch beyond the transition point and have to be replaced by the coexistence curve, e.g. via the Maxwell 
construction. Along these branches the metastable state is thermodynamically stable as long as the second free 
energy derivative is positive and the escape from from the local minimum takes places via thermal activation. The 
theory of Ginzburg and Landau takes into account localized fluctuations, but these cannot be very large, excluding 
therefore interesting regions from its range of validity, like the one around the critical temperature. At the end of the 
metastable branch is the spinodal point. There the second derivative of the free energy changes of sign, making the 
whole system unstable. The decay to the stable phase is global and quick. Beyond mean field, statistical mechanics of 
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short-interacting systems cannot account properly for metastability, since the partition function in ensemble theory 
is dominated by the global minimum of the free energy functional in phase space. 

Metastable states with infinite lifetime in short range system with translational invariance are thermodynamically 
forbidden by stability requirements [28| , although they can be obtained by setting specific constraints to the accessible 
phase space or taking suitable limits in some perturbative approaches The constraint is related to a timescale, up 
to which the metastable phase do not explore the whole available phase space. In this context, in the droplet theory 
for the condensation of the field-driven Ising transition Q , the fluctuations neglected in mean-field approaches give 
raise locally to droplets of the stable phase. An ensemble of non-interacting droplets of the stable phase is considered, 
its growth being favored by the bulk free energy, which is lower in presence of an external field, and hindered by the 
surface tension. The model has been extensively tested in two and three dimensions Q- 

For the T-driven Potts model transition, differently from what happens in the Ising case, the growth of droplets 
just below T q is favoured by the tendency to lower the total energy of the system (the interface perimeter between 
droplets), to the detriment of the system entropy (which is higher when many small droplets are present). In this 
case, the ensemble of droplets is also different with respect to the Ising case in that droplets have an intrinsic entropy, 
and can assume different shapes. Meunier and Morel [ijj] considered an adaptation of the droplet approach to the 2D 
Potts model. Their theory leads to a temperature range in which metastability is present in finite-size systems, range 
which, however, shrinks to zero in the thermodynamic limit. This behaviour, related in some way to the anomalous 
growth of fluctuations with the system size near the critical point, has, to our knowledge, not yet been confirmed 
by dynamic calculations. Investigation of metastable states from a pure dynamical point of view seems therefore to 
deserve interest. 



Numerical simulation of the Potts model with more than four states in 2D constitutes a very stimulating and 
effective way for investigating the nature of metastable states in systems with many low-temperature stable phases. 

A way for identifying metastable states just below the transition temperature is to look at the time behaviour 
of the energy per site e(t). As a test system we have considered the Potts model on square lattices with q = 
12, 50. Simulations have been performed by preparing the system in the high temperature disordered phase and then 
quenching it by Metropolis dynamics below the transition temperature T q . 

For the q = 12 system, for which = T q = l//3 q = 0.668, we have taken the quench inverse temperature at 
(3 = (3 q (l + j ■ 0.005), with j = 1, . . . 10. The resulting time behaviour for lattices of size L = 250 is shown in figureQ] 
Each line is obtained by averaging over 50 different realizations of initial conditions and thermal noise, and the errors 
are the resulting standard deviations. It is seen that for temperature close to T qi e{t) stays stationary (a signature of 
metastability [J]) until large times, at an excess energy with respect to the equilibrium state. The existence of these 
metastable states is in agreement with [lfjj . who however predict their disappearance in the thermodynamic limit. 

From the point of view of the local dynamics, if the system is in a metastable state just below the transition 
temperature, it is expected to be indistinguishable from one in equilibrium. This property is reflected in the time 
behaviour of the self-correlation function. In the equilibrium disordered phase (T > T q ) this function displays a fast 
(exponential) relaxation. The same behaviour is therefore expected for the metastable states below T q , but will no 
longer be true after the decay towards the stable state has begun. The correlation function computed waiting different 
times t w after the quench time should therefore signal this decay through the onset of a slower decay (ageing) . The 
self-correlation function can be computed as: 



where t w is the time elapsed after the quench before starting the computation and the sum is over all the lattice site 
(averages over different realizations of initial conditions and thermal noise are taken). 

As mentioned before, it is expected that for t w short enough the relaxation is fast. Figure [2] shows c(t,t w ) for 
different t w for a system with q = 50, L = 100 and T = 0.470, being T q = 0.479. Curves are obtained by averaging 
over 10 realizations (except one indicated in the legend). It is seen that after a very short transient time c(t,t w ) 
becomes translationally invariant for long times, also showing no dependence on the size in the investigated range. 

Figure [3] shows c(t,t w ) vs time t for q = 50, t w = 1000 and different temperatures. Ageing becomes larger when 
the temperature is decreased, showing that the system is already relaxing toward the ordering phase. 

These results show the existence of well definite metastable states on square lattices of finite size. The question 
of their existence in the thermodynamic limit is thus open. Meunier and Morel fl(i| have recently faced the problem 
by investigating the energy probability distribution at the transition temperature, concluding that metastability is 
observable only in systems of finite size. The anomalous growth of energy fluctuations near the critical point is at 
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the origin of this fact, but the microscopic mechanisms behind this shrinking of the metastability interval for large 
systems remain unknown. It would be very interesting, however, to have a description of this phenomena in terms of 
microscopic clusters. One could then compare the mechanisms of nucleation in the Potts case with those predicted for 
the field-driven Ising case by Classical Nuclation Theory. In Classical nucleation theory, when the system size is much 
larger than the lenthscales involved in the the nucleation process then lifetimes and properties of the metastable phase 
do not depend on the system size Q. What is, hence, the microscopic difference in the T-driven Potts transition, 
with respect to this scenario, which leads to the shrinking of the metastable temperature range for large system? 

According to the nucleation theory the finite lifetime is usually determined by the time needed to the system 
for nucleating a drop of the stable phase bigger than the critical size. At the same time, as mentioned above, the 
metastable state displays short relaxation times similar to the high temperature phase. The latter usually increase 
when temperature decreases, and it may happens to become equal to the nucleation time at a certain temperature 
T s . This temperature corresponds to a pseudo-spinodal point [29( and sets a limit to the observability of undercooled 
metastable states. The first point is thus to establish if T s < T q even for infinite systems or if the two temperature 
merge. Another point concerns the growing dynamics following the nucleation. In the Ising model the metastable 
state nucleates a critical droplet that afterward grows very fast because of its favourable free energy with respect 
to the unstable phase. In the temperature driven Potts model however each phase is equivalent to the others and 
competes with them for the growth. 

IV. CONCLUSIONS 

We have discussed some aspects of the off-equilibrium dynamics of the q > 4 Potts model quenched below the 
transition temperature. We believe that this can be relevant for understanding the non-equilibrium and the ordering 
dynamics of systems below a first order phase transition temperature, in the presence of many competing equivalent 
ground states. After recalling some recent findings we have focussed on some characterization of the undercooled 
metastable states which are observable after a quench below the transition temperature. These states display a finite 
lifetime which depends on the temperature and are detectable through the existence of a plateau in the energy (well 
above the equilibrium energy) and the fast relaxation time. Being their existence in the thermodynamical limit 
questioned by recent work, more extensive numerical simulations will be necessary to support this statement. More 
work in this field is highly desirable also for understanding the detailed mechanism through which a system starts to 
order after the quench in the presence of many degenerate states, which is not yet well established. 
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FIG. 1: Time behaviour of the energy per site, after the quench of the model with q = 12 at different quenching temperatures 
below T q (see text). The quenching temperature increases from bottom to top. 
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FIG. 2: Time behaviour of the self-correlation function for different waiting times for a system with q — 50. The lattice is 
L = 100 in size and the quench is at the temperature T = 0.470 (T q = 0.479) 
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FIG. 3: Time behaviour of the self-correlation-function for t w — 1000 after quenching at different temperatures below T q = 0.479. 
At low temperatures ageing appears. 



